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Abstract. We use quantum graphs as a model to study various mathematical aspects of 
the vacuum energy, such as convergence of periodic path expansions, consistency among 
different methods (trace formulae versus method of images) and the possible connection 
with the underlying classical dynamics. In our study we derive an expansion for the vacuum 
energy in terms of periodic paths on the graph and prove its convergence and smooth 
dependence on the bond lengths of the graph. For an important special case of graphs 
with equal bond lengths, we derive a simpler explicit formula. With minor changes this 
formula also applies to graphs with rational (up to a common factor) bond lengths. The 
main results are derived using the trace formula. We also discuss an alternative approach 
using the method of images and prove that the results are consistent. This may have 
important consequences for other systems, since the method of images, unlike the trace 
formula, includes a sum over special "bounce paths" . We succeed in showing that in 
our model bounce paths do not contribute to the vacuum energy. Finally, we discuss the 
proposed possible link between the magnitude of the vacuum energy and the type (chaotic 
vs. integrable) of the underlying classical dynamics. Within a random matrix model we 
calculate the variance of the vacuum energy over several ensembles and find evidence that 
the level repulsion leads to suppression of the vacuum energy. 
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1. Introduction 



Vacuum energy is a concept arising in quantum field theory and was first shown by Casimir 
[1] to have an observable efi^ect on two perfectly conducting parallel plates, causing them to 
attract. Since then, experiments with various physical geometries have confirmed the effects 
of vacuum energy (see [2, 3, 4, 5]). 

In time-independent situations the vacuum energy is formally given by 



where k'^ are the eigenvalues of a Hamiltonian, H, where throughout we take h— 1 — c. The 
above expression arises in quantum field theory in the context of cavities and cosmological 
models [3], and it is formally divergent. To get meaningful result from this expression, the 
vacuum energies for two different configurations are subtracted from one another [2]. To 
accomplish this in a systematic way, we employ an ultra-violet cutoff defining the energy as 
the regular part of 



as i — > 0. To evaluate (2), it is sometimes convenient to employ the trace of the cylinder 
kernel, T{t) = E„e~''"*, [6]. In this way, E{t) = —T'(t)/2. The singular term in the 
expansion of E{t) is related to the vacuum energy density of free space, and physical 
justification for its removal is described, for example, in [2] (for systems similar to those 
considered here, see [7, 8]). 

A widely employed method of calculation of the vacuum energy is expanding it into a 
sum over classical paths [9, 10, 11, 12, 13, 14]. The expansion is usually done by the method 
of images, or "multiple refiections", leading to a sum over all closed paths. It has been 
argued in [11] that for certain geometries restricting the sum to include only the periodic 
paths ( "semiclassical evaluation") correctly reproduces asymptotic behavior of the vacuum 
energy and is much simpler to evaluate. A periodic path comes back to the starting point 
with the same momentum, while a closed path might not. Another popular approximation 
predicts the sign of the vacuum energy by considering only short orbits [15, 8]. This implicitly 
assumes that the convergence of the complete series is sufficiently fast. 

In the present paper we aim to contribute to this discussion by studying the vacuum 
energy on quantum graphs (for another model where similar questions are addressed, see 
[16]) . Quantum graphs are often used as mathematical models that exhibit the relevant 
phenomena while being sufficiently simple to allow mathematical treatment. We compare 
the method of images with the direct application of the trace formula (which is exact on 
graphs) and demonstrate that the outcome is the same. This is done by showing that the 
contribution of the "bounce paths" — the paths that are closed but not periodic — is 
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identically zero. We also prove that the resulting sums converge, giving an estimate for the 
rate of convergence, and can be differentiated term-by-term with respect to the topological 
parameters present in the model. 

One of the main reasons for the success of quantum graphs as models (for example, of 
quantum chaos, see [17] for a review) is the existence of an exact trace formula on quantum 
graphs. A trace formula is a relation between the spectrum and the set of periodic orbits of 
the system. For graphs, the trace formula was first found by Roth [18] and then by Kottos 
and Smilansky [19]. Subsequent studies of the mathematical properties of the trace formula 
on graphs included works by Kostrykin, Potthoff and Schrader [20] and Winn [21]. 

The trace formula of [19] gives an expression for the density of states d{k) defined as 

oo 

d{k) = Y,6{k-k^), (3) 

n=l 

where S{-) is the Dirac delta-function. The vacuum energy is then, formally, 

E{t) = ^ y ke-^^d{k)dk, (4) 

quickly leading to the result (for fc-independent scattering matrices, see Sec. 2), 

1 ~ A 

n=l peV P P 

The sum is over the set of all periodic paths V on the graph, £p is the metric length of the 
path p, Up is the period (the number of bonds) of the path and Ap is its stability amplitude 
(see Section 4.1 for definitions). 

After introducing the notation in Section 2 and considering some explicit examples in 
Section 3, we prove the mathematical correctness of the calculation outlined above. This is 
done in Section 4.3, where the convergence of (5) is also analyzed. In Section 4.4 we show 
that we can differentiate (5) with respect to individual bond lengths, showing the smoothness 
(C°°) of Ec as a function of lengths. 

In Section 5 we briefiy discuss the method of images (covered more fully in [22]) and 
show that the contributions from the bounce paths cancel. Finally, in Section 6.1 wc discuss 
random matrix models of the vacuum energy. As expected, in such models the average 
vacuum energy is zero and there is no preferred sign to the Casimir force. However, by 
analyzing the second moment of the energy we confirm an earlier observation by Fulling 
[23, 6] that level repulsion tends to decrease the magnitude of the energy. 

2. Vacuum energy and quantum graphs 

Quantum graphs were introduced as a model of vacuum energy by Fulling [24] who considered 
the effect of the energy density near a quantum graph vertex by constructing the cylinder 
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kernel for an infinite star graph (a graph with one vertex and B bonds extending to infinity). 
The quantum field theory origins of this in a graph context were given by Bellazzini and 
Mintchev [25]. The vacuum energy expression for quantum graphs obtained in the present 
manuscript was also used in [8] where the convergence was investigated numerically (we 
prove rigorous estimates here). 

We start by briefiy recalling the terminology of the quantum graph model, see [26] for 
a general review of quantum graphs. We consider a finite metric graph F consisting of a set 
of vertices V, and a set of bonds B. A (undirected) bond b connecting the vertices v and w 
is denoted by {v,w}. Each bond b is associated with a closed interval [0, L;,], thus fixing a 
preferred direction along the bond (from to L^) . This direction can be chosen arbitrarily. 
If the direction from f to w is chosen, the bond b = {v,w} gives rise to two directed bonds, 
= {v,w) and = (w.v). Whenever the distinction between 6+ and b^ is unimportant, 
we will denote the directed bonds by Greek letters: a, j3. In addition, the reversal a is 
denoted by a (e.g. if a = then a = b~). We will denote by B the number of bonds \B\; 
correspondingly, the number of directed bonds is 2B. The length of the directed bond is 
naturally determined by the length Li, of the underlying undirected bond. The total length 
of r is £ = SbGS ^b- We will denote by L = diag{Li, . . . , Lb, Li, ... , Lb} the diagonal 
2B X 2B matrix of directed bond lengths. 

In this article we study the spectrum of the negative Laplacian on the graph. The 
Laplacian acts on the Hilbert space 7i(r) := ^^g^j -^^^([0, Li,]) of (Sobolev) functions defined 
on the bonds of the graph. On the bond b it acts as the 1-dimensional differential operator 
— A domain on which the Laplacian is self-adjoint may be defined by specifying matching 
conditions at the vertices of F, see e.g. [27, 28, 29, 26]. 

To specify the matching conditions, let / be a function in 7i(r). For a vertex v of 
degree d we denote by f^"^ the vector of values of / at v, f^''^ = {fbi{v), . . . , fbd{v))^, where 
fb{v) = /fo(0) if 6 = {v, w} is oriented from vtow and fb{v) = fb{Lb) otherwise. Furthermore, 
let g^^) denote the vector of outgoing derivatives of / at v, g^'"^ — {fl^{v), . . . , fl^{v)Y , i.e. 
fl{v) — fl,{0) lib — {v, w} is oriented from v to w and fl{v) — —fl,{Lb) otherwise. Matching 
conditions at v can be specified by a pair of matrices A^") and B^^^ through the hnear equation 

A('')fM +B('')g('') = . (6) 

The matching conditions define a self-adjoint operator if (A'^'^B*^*)) has maximal rank and 
A^'^B*^*)''' is self-adjoint at each vertex (where a B*^*)''' represents the adjoint of B). 
A solution to the eigenvalue equation on the bond b, 

d^ 

-^'06(a;6) = k'^ipb{xb), (7) 

can be written as a linear combination of plane waves, 

^b{xb) = Cbe'^''' + Cbe-^^''' . (8) 
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where c is the coefficient of an outgoing plane wave at and c the coefficient of the incoming 
plane wave at 0. A solution on the whole graph can be defined by specifying the corresponding 
vector of coefficients c = (ci, . . . , cb, ci, . . . , ■ 

The matching conditions at the vertex v define a vertex scattering matrix 



(k) = -(A^'') + im(''))-^(A(T;) - iA;B(")), 



(9) 



see [27]. a'^'"^ is unitary and the elements of cr^'"^ are complex transition amphtudes which 
in general depend on k. However, for a large class of matching conditions including the so 
called Kirchhoff or natural conditions the S-matrix is independent of k. Kirchhoff matching 
conditions require that ip is continuous at the vertex and the outgoing derivatives of ip at 
the vertex sum to zero. These conditions may be written in the form (6) with matrices 
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(10) 



Substituting in (9) leads to /c-independent transition amplitudes 

2 

Mij = ^ - ^ij , (11) 

where d is the degree of v. 

The matrix cr^^^ relates incoming and outgoing plane wave coefficients at v, c^""^ = ac^'"\ 
Collecting together transition amplitudes from all the vertices of a graph we may define the 
familiar 2B x 2B bond scattering matrix S [30], 



[^]{v',w'){v,w) — Sw,v'[<T^'^'](v^w') ■ (12) 

Wc shall also need the quantum evolution operator U = Sc^'^'", which acts on the vector of 
2B plane wave coefficients indexed by directed bonds. For a general graph, the spectrum 
can be computed as the zeros of the equation 

det(I - Se*^^) = 0. (13) 

This formula goes back at least to [31]; for a discussion of scattering matrices of different 
types we refer the reader to [30]. 

The spectral theory of quantum graphs is often extended to included quantum evolution 
operators defined by specifying a priori a set of unitary vertex scattering matrices as in 
[32, 33]. The vertex scattering matrices are typically chosen to be k independent and to 
have other desirable features, for instance transition amplitudes of equal magnitude as is the 
case if the scattering matrix is a fast Fourier transform matrix 



1 12:^17 

— e d 



Vd 



(14) 
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Such a scattering matrix would be difficult to produce from matching conditions of a self- 
adjoint operator at a vertex. For such a quantum evolution operator the spectrum is still 
defined by (13) however the scattering matrices in general no longer correspond to a self- 
adjoint realization of the Laplace operator on the graph. 

In the following we work with fc-independent scattering matrices which can either be 
considered to come from a self-adjoint Laplace operator or through specifying unitary vertex 
scattering matrices directly. 

3. Some explicit examples 

3.1. Star graph with bonds of equal length 

One case where the vacuum energy can be computed explicitly is the quantum star graph with 
bonds of equal length. This example was first considered in [8] and here, for completeness, 
we summarize the computation. 



Consider a star graph (see Fig. 1) with B bonds. The bond h has length Li,. We consider 
the equation (7) with Neumann conditions. At the central vertex, this translates into 




Figure 1. A star graph with B bonds 
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(15) 
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(16) 
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Z{k) = ^tan(A;L6) = 0, 



6=1 
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if the lengths are rationally independent. We note that since Z{k) is an increasing function, 
there is exactly one zero of Z{k) between each pair of consecutive poles. Because of rational 
independence, the poles of different tangents do not coincide. If we lift the restriction on 
lengths, in addition to the zeros of Z{k) we have the following eigenvalues: k, is an eigenvalue 
of multiplicity m if there are m + 1 lengths L^^ such that k is a pole of each tan(A;Lf,^). 

In particular, if all lengths are equal, Li = . . . = Lb = L, Z{k) is simply Btaii{kL). 
Each zero of Z{k) is a simple eigenvalue, while each pole is an eigenvalue of multiplicity 
B — 1. We can now evaluate 

T{t) = e-*'" = + - 1) I] e-(2-+i)-*/2i ^ ^_ -it/L 

n=l n=0 ?i=0 

BL . 1 (B-3)7r ^,2, 

Now we take the regular part and the limit t — > of —T'{t)/2 which gives us the vacuum 
energy, 

Bc = ^^^. (17) 

48L ^ ' 

3.2. General graphs with bonds of equal length 

If all bond lengths of the graph are equal to L, we can use equation (13) to explicitly 
describe the infinite spectrum of the graph in terms of the finite spectrum of S. Indeed 
let A be the diagonal matrix of the eigenvalues of S. Then e*^^ = e^'^'^I and, therefore, 
det(I — Se*'^''"^) = det(I — Ae'*^^). Consequently the solutions of equation (13) are the values 
k such that 

for some j. Here by e'^^ we denoted the j-th eigenvalue of (unitary) matrix S. Thus the 
/c-spectrum is 



n=l 



where we choose 6j to lie between and 27r. 



Now we compute the trace of the cylinder kernel T(t) 

oo 2_B oo 2B 



n=l ■j=\ n=\ j=l 

" L e.-TT 3^2 + 27r2 - 66^71 



= E 



2Tit 2ti 12Lti ^ ' 



(21) 
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Thus, the vacuum energy is 

Ec^-Yl ' i2L7r ^ "L ^ Md^M, (22) 

i=i i=i 

where B2{-) is the second Bernoulh polynomial (compare to the B = 1 case discussed in [6]). 

As an example of application of formula (22) consider again the star graph with equal 
bond lengths. The eigenphases 9j of the S-matrix of a star graph are 0, tt, it/2 and 37r/2, 
the latter two with multiplicity B — 1. Substituting into equation (22) one can recover (17). 

3.3. Graphs with bonds of rational length 

By introducing the Neumann vertices of degree 2 we do not change the spectrum of the 
graph and therefore the vacuum energy. On the other hand, if the bonds of the graph are 
rational (up to an overall factor), by introducing such "dummy" vertices we can convert 
the original graph into a graph with bonds of equal length. The number of bonds (and the 
dimension of the scattering matrix S) will increase as a result, but the vacuum energy will 
still be explicitly computable using equation (22). 

Moreover, one can conceivably approximate rationally independent lengths by rational 
ones and use the result as a numerical approximation to the true vacuum energy. For this 
approach to work one needs to know, a priori, that Ec is continuous as a function of bond 
lengths. This question is one of the main subjects of Section 4. 

4. Vacuum energy via the trace formula 
4.1. Formal calculation 

In this section we perform a formal calculation of the vacuum energy Ec using the trace 
formula. We shall investigate the rigor of the manipulations in Section 4.3. 

The trace formula (see, e.g. [17]) connects the spectrum {kn} of the graph with the set 
of all periodic orbits (or periodic paths) on the graphs. A periodic path of period n is a 
sequence (cti, «2, . . . , «„) of directed bonds which satisfy [SJc^^-^^a^. 7^ for all j = 1, . . . ,n 
(the index j + 1 is taken modulo n). A periodic orbit is an equivalence class of periodic paths 
with respect to the cyclic shift (ai, 0:2, . . . , ctn) ^ {c(2, • • • , «n, cti). We denote by Vn the set 
of all periodic paths of period n and by V the set of periodic paths of all periods. The the 
trace formula can be written as 

°° r 1 f ■ 

d{k) = J2S{k-kn) = - + -Re J2 ^p— e'''^^ (23) 

n=l ^ ^ pGV 

where £ is the total length of the graph (the sum of the bond lengths), Up is the period 
of the periodic path p, ip = Y^jLi^aj is the length of p and Ap = YYj=i[^]aj+i,aj is its 
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amplitude. Using the trace formula and equation (2) we can formally compute the vacuum 
energy. Indeed, 



oo 



VA;„e-**="=/ ke-^'d{k)dk (24) 

n=l 



27r Jo 



y ke-'''dk + -ReY,\^ J ke-'^'+'^'^-dk (25) 



per -'f-^o 



Removing the divergent Weyl term C/nt^ due to regularization and taking the limit t — >■ 
leads to the following simple expression for the vacuum energy, 

E,^-^Rey^. (27) 
4-2. Equal bond lengths; equivalence to (22) 

In the case of equal bond lengths (27) should be equivalent to the sum of second Bernoulli 
polynomials (22). If the length of each bond is L an orbit that visits n bonds has length nL 
and we may rewrite (27) as a sum over the topological length n followed by a sum over the 
set of all periodic paths visiting n bonds, Vm 

1 °° 1 

^^"=-^R«E^E-4. (28) 

n=l p&Vn 

^ OO ^ 2B 2B 

— ~ 2tiL n? ^ ^ ' ' ' ^ ^ ^^{^aia2^a2a3 ■ ■ ■ ^a„ai) (29) 

n=l oi=l Q;n=l 

1 °° 1 

2.TE^(t^S" + U(S<r) (30) 



n=l 

oo 2B 



_ 1 COS nOj 

n=l j=l 

where e'^^ are the eigenvalues of the matrix S with ^ ^ ^ 27r. The sum over n can be 
expressed in a closed form, see Abramowitz and Stegun [34], formulae 27.8.6, 

V = ^0^ + 2.^ -^^0 ^ (32) 

n=l 

where i^2(-) is the second Bernoulli polynomial. Consequently we recover expression (22). 
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4-3. Convergence of (27); vacuum energy as a function of the bond lengths 
We shall now present a rigorous derivation of equation (27). 
Theorem 1. The vacuum energy of the graph, defined as 



Ec^- lim V knc-*^"^ - — , 

n=l 



is given by 



1 °° 1 poo 



(33) 




pflp 



(34) 



where Vn denotes the set of all periodic paths of period n. The vacuum energy is smooth 
(C°°) as a function of bond lengths on the set {L^ > 0}. 

Remark 1. The sum over the periodic orbits in (34) is finite for each n. We will show, in 
particular, that the sum over n is absolutely and uniformly convergent. More precisely we 
will derive the following bound. 



where 2B is the number of (directed) bonds and Lmin is the minimal bond length. This 
estimate shows that, if the (finite!) sum over periodic orbits of a fixed length is performed 
first, the series in (34) becomes absolutely convergent. Moreover, it is uniformly convergent 
with respect to the change in bond length as long as L„^in remains bounded away from zero. 

We would like to mention that our estimate (35) agrees with the numerical results of [8], 
even though in [8] the ordering of the periodic orbits was different (according to the metric 
length ip rather than topological length n) . 

Proof of (33)-(34)- The C°° part of the proof will be given in the following section. 
We start with the definition of the vacuum energy and integrate by parts. 



where N{k) is the integrated density of states (IDS), apiecewise constant, increasing function 




(35) 




(36) 



N{k) = #{n : < A;„ < A;}. 
The integrated density of states N{k) can be spht into two parts, 

N{k) = const + — + N'^'^ik). 



(37) 



(38) 
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The first two terms are unimportant: The first term makes no contribution in the integral, 
and the second term is removed at the regularization stage. The oscillatory part possesses 
an expansion, see [17], equation (5.24), 



N° 



1 °° 1 

(A; + i£) = -Im^-trU'^(A; + i£), (39) 



n=l 



where U = Se^^^^. 

This expansion is absolutely convergent as long as £ > since the matrix U"(A; + ie) is 
then subunitary (all eigenvalues he within a circle of radius strictly less than 1). As £ ^ 0, 
N'^'^{k + ie) converges to N°^^{k) pointwise almost everywhere. Moreover, \N°^^{k + ie)\ is 
uniformly bounded by the number of bonds B (in other systems one can show that the Weyl 
law implies that \N°^^{k)\ — 0{k'^) as /c — > oo where d is the dimension of the system). 
Therefore, 

1 

E^^-- lim lim / {tk - l)e-*''N°'%k + ie)dk, (40) 
and, using the convergence of expansion (39), 

1 ^ 1 roc 

-— limlim V-Im / (tfc - 1)6"*'^ [tr U"(A; + ie)] dA;. (41) 

n=l " 

We will now show that the integral 

POD 

Rn= {tk-l)e~'^[tvV''{k + i£)]dk (42) 
Jo 

is absolutely bounded by 1/n. Thus the series is absolutely convergent uniformly in e and t 
and we can take the limits inside the sum. 

A typical term in the (finite!) expansion of the trace is Ape^''^''e~^^p . The two exponential 
factors, e"*'^ and e^'^^'', ensure that the integrand is exponentially decaying in k in the first 
quadrant of C. Therefore we can rotate the contour of integration to the imaginary line, 
k — is. The integral becomes 

poo 

Rn^i {ist - l)e-i^* tr (Se-(*+^)^) " ds. (43) 
Jo 

We estimate 

2B 

tr (Se-(*+'^)^)"| < 5^ |A,-(s + e)|", (44) 
i=i 

where 2B is the size of the matrix S and \j is j-th eigenvalue of the matrix Se^*^*^^-*^. 
According to a familiar argument (see, e.g. [35]), the maximal |Aj| is bounded from above 
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by the maximal singular value of the matrix. The singular values are square roots of the 
eigenvalues of 

which is a diagonal matrix with the maximal entry e"^*^*"*"^^^™'" (-^>min is the smallest bond 
length of the graph) . Thus we can estimate 

tr(Se-('+")^)" <e-""^-'°2B. (45) 

Finally, 

\RJ < 2B / list - lle-^^-^-'-c^s . (46) 

Jo n 

We notice that the integrand of (43) can be absolutely bounded by e~"^'^*^'^™'°~'') for 
arbitrarily small 5 and sufficiently small t. Thus, having brought the limits inside the sum, 
we can use the dominated convergence theorem to bring them inside the integral. Taking 
the limit t — > and £ — > inside integral (43) produces 



1 °° 1 /"OO 

E.^—T-Rc tr(Se-^^)"ds. 

2^ n=l ^ 



(47) 



Now we expand the trace, 

tr (Se-^^)" = Yl ^pe-^^^ (48) 

and integrate term by term to recover (34). □ 

Remcirk 2. The basic idea of the proof, shifting the convergence into the subunitary matrix 
(see equation (43) and (45)) can also be used to study the convergence of the trace formula 
itself. This was done in [36]. 

RemcLrk 3. Since S is /c-independent, the trace in the integral is real and we do not need 
to take the real part in equations (33)- (34). Indeed, the following lemma easily follows from 
the definition of S and [20, Prop 2.4] appfied to the matrices cr^'"\ 

Lemma 1. The S-matrix of a graph is k-independent if and only if it satisfies 

jsj = s^ (49) 

where J is defined by Ja,/3 — S^^^ (P is the reversal of j3 as defined in section 2). 

Now the complex conjugate of tr (Se~*^)" is 

tr (e-^^St)" = tr (e-^^JSJ)" = tr (SJe-^^J)" = tr (Se"^^)" , (50) 

where we used the fact that the length of a bond is invariant with respect to direction reversal 
and, therefore, Je~*^J = e~*^. 
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4.4- Derivatives of the vacuum energy 



Proof of dijferentiability of Ec- We differentiate the expression (33) term by term and show 
that the result is also absolutely convergent. This can be done by using the following bound 
on the partial derivatives of tr (Se"*^)", 



tr (Se- 



(i^min/2)l™l ■ 



Before proving (51) we note that it implies the following bound, 



Consequently 



1 

n .10 

OO 

E- 

^ n 



tr (Se-^^)''ds 



9L^i . . . 9L^^ 



2B 



n2(W2)IH+i ■ 



(51) 



(52) 



Qmi QraB 



converges absolutely and we arc therefore allowed to differentiate (33) term by term. We 
conclude that the vacuum energy is as a function of bond lengths. 

To prove bound (51) we use the Cauchy integral formula (see, e.g., [37]), 



tr (Se-^)"ds 



Qfmi QrriB 



sL\n 



■tr(Se-"^) 



1 



(2vr) 



B 



•2w 



2- tr(Se-^(L+^('^)))" 



.7 



d(f)i . . . d(f)B 



(53) 



B 



where Rj = r^e^^^ and R(0) = diag{i?i, . . . , Rb, Ri, ■ ■ ■ , Rb}- Let A = Sc-*(^+^). Since S 
is unitary, we find A^A = e^''(^^+^+^^). The eigenvalues of A are bounded by the maximal 
singular value of A, 

|eig(A)| ^ max e-^(^''+'-''=°'^('^'')) . (54) 

6=1.. .B 

By choosing the radius rf, — L^am/'^ we see that 



|eig(A)| ^ max e 
' ^ ^' 6=1. ..s 



./2) 



and 



2B 



|tr(Se-"(^+^('^»)"| = |tr A"| < |eig(A)jr < 2Be~'''^^'-'^ . 
Using this bound in (53) gives 



gmi QmB 



tr(Se-^^)" 



9L^i . . . 9L^^ 
which establishes (51). 



(2 



1 r 
WL "Jo 



(i^min/2)l™l 



d0i . . . d4>B 



(55) 
(56) 

(57) 
□ 
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5. Method of images expansion and the equivalence of two expansions 



We can evaluate the trace of the cyhnder kernel T(t) by constructing the kernel itself. The 
cylinder kernel, Tw{t]x,y) where x is measured on bond h and y on bond h' satisfies the 
following equation on each bond h 

92 , ^ 



Tw {t; X, y) = -^T^ {t; x, y) (58) 

for i > 0, Tftft/ — > as t — > 00, with boundary conditions (6) with respect to the x coordinate, 
and the initial condition Tw{Q] x, y) — 5w^{x — y)- By separating variables in (58) it can be 
shown that the trace of T is 

T{t) V / {t; x,x)dx = T e-'^"* . (59) 
6=1 

The corresponding free space kernel (i.e. the solution of equation (58) on the whole real 
line with the initial condition To(0; x — y) — S{x — y)) is 

We can apply the method of images (i.e. multiple reflection) to the free space kernel to obtain 
the kernel on the graph. The cylinder kernel is then written in terms of paths which begin 
at y on the bond h' and end at x on the bond h (the details of the construction are given in 
[22]), 

Tw{t\ X, y) = hyToit] x - y) 

00 

+ X X [^b+pb'+^o(^; L^, + ^p + x-y)+ ^b+p;,'-To(t; -^p-x- y) (61) 

n=0 p€P„ 

+^h-pb'-To{t] -ip - Lh + X - y) + Ah-pb'+To{t; Lh> + £p + U - x - y) 

In the above expression, we denote the two directed bonds associated with the undirected 
bond h by and h~ . The path of topological length n, p = (ai, . . . , a„) is an n vector 
of directed bonds and P„ is the set of all such paths. The metric length of a path 
'P = YTi=iLaj and Ap = [S] ■■• [Sjagaa [S]„2«i is the stability amplitude of the 



IS 



path. Note that even at the stage represented in (62) we have assumed the matrix S is 
A;- independent. 

To obtain the trace, we let b — b' and y — x. While this corresponds to "closing" the 
paths, we do not always get periodic paths in the topological sense of Section 4.1. Indeed, a 
periodic path would return to the initial point Xf, with the same momentum (or direction), 
whereas when the paths corresponding to the second and fourth terms of (62) return to Xb, 
the momentum has the opposite sign, see Fig. 2. The latter paths we shall call bounce paths. 
The difficulty in the method of images is in the handling of the bounce paths. 
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Figure 2. An example of a periodic path (left) and a bounce path (right). 



After integrating Ti,h{t; x, x) we break the formula into three parts 

Tit) = Tps(t) + Tpo(t) + TbpW, (62) 

where FS stands for free space, BP for bounce paths, and PO for periodic orbits. The three 
parts are (taking into account that To(t;a;) is even in x): 

^FS W = E / 0) = 0)^' (63) 



=1 

B 







f r 1 

^POW = H HY. [A+p6+ro(t;i>b + ^p) +^-p6-To(t;£p + Lb)J dx, 

6=1 "^0 n=0 peP„ 



B „L oc 



(64) 



Y [A+pb-To{t;ip + 2x) + A-p6+ro(i; + £p - 2x)J dx. (65) 

We shall use the following lemma to simplify the bounce path term. 
Lemma 2. // the scattering matrix S of a graph is k -independent, then 

2B 
a=l 

where Ja,(3 = 

Proof. Writing out the above, 

^ ^ ^aa„-i---aia — 'S'a„_ia„_2 ' ' ' Sa3a2^a2ai ^ ^ SaiaSaa„-i- (66) 

The last sum in the above is equivalent to an element of SJS, but JSJS = I, see lemma 1, 
and J-^ = J, therefore SJS = J. □ 

Now we come to the following theorem which proves the equivalence of this method to 
the trace formula method of Section 4. 
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Theorem 2. 

T(^) = ^ + itr(J5) + f:$:A,§^-^ (67) 



n=l peVn P 

and consequently, 

Ec^-^yy ^. (68) 

n=l pGVn ^ 

Remark 4. A special case of this proof was done for the heat kernel with Kirchhoff conditions 
by Roth [18]. Similarly, Kostrykin and Schrader have an analogue of this proof (again, for 
the heat kernel) in 



Proof. The first term in T{t) is the free space term Tpg, which we found above to be equal 
to To(t; 0)C = C/nt. Thus we need only consider the periodic path and bounce path terms. 
Periodic orbit contribution. There is no a;-dependence in (64), so the integration gives 

oo B 

Tpoit) = E E E K-P"- + A+pb+]To{t]ip + L,)L,. (69) 

n=0 peP„ b=l 

The two amplitudes can be written as a single sum over the directed bond a, 

B 2B 

J2 [^&-P&- + A+P6+] ^o(^; £p + LbjLh = A^pMt; + La)L^ . (70) 

6=1 a=l 

The path apa is actually a periodic path of period n + 1 which we will denote by p. The 
amplitude of p is Ap = A^pa and its length is ip = ip + La- Thus, 

oo . oo . 

Tpoit) -EE Ap^T,{t;ip) = E E \iTo{t;ip) . (71) 

n=0 peV„+i n=l peVn 

Bounce path contribution. Wc can again combine the two amplitudes and change the 
variables in the integrals to obtain 

^BpW = 2EEE^"pW Ut;x)dx. (72) 

n=0 peP„ a=l "'^P 

We now fix 1/ > and introduce the cutoff function, 

H{y-x) = { . (73) 

10 X > y 

Instead of To(t; x) in the integral, we consider Ty{t] x) ~ To[t; x)H{y—x). Since the minimum 
bond length Lmin is greater than zero, taking m large enough we will have ip > y for all 
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paths of topological length m — 1 or greater. Therefore, for a path p in P,„_i or in and 
any a, we can write 



/ 



Ty{t;x) dx = Ty{t;x) dx, 



(74) 



since the integrand is identically zero on both intervals of integration. We can also ignore 
all paths from P„ with n > m. 

We can therefore write the bounce path contribution with the given cutoff in the 
following form, 

Ty{t; x) dx 



n=0 peP„ a=l 

+ 2 ^ 1]^«P" / fy{t]x)dx^- J] J^^apa / fy{t]x)dx 



(75) 



Applying Lemma 2 to the last set of sums in (75), we obtain 



2B 



2S 



Tj,(t; x) dx. 



(76) 



Tj,(i;x)dx= J^^/3p/3 
peP^_2 /3=i 

Here the new path is the same as the old path but with the first and the last bond removed 
and 13 corresponds to the last bond of the old path. 

Now we can take the sum corresponding to n = m — 2 in (75) and add it to the result 
of (76), 

2S .£p+2L„ 2B 



pGPm-2 a=l 

2B 



fy{t-x)dx+ 

peP^-2 /3=i 



.+2Lf} 



Ty{t] x) dx 



= E E^«P« fTy{t;x)dx (77) 

p€Pm-2 a=l ^P 

Therefore, Tgp(t) can be rewritten exactly in the form of (75) but with m reduced by 1. 
Proceeding by induction, we obtain 

a "^O q;=1 "^^f 

However, J^q = and we can take the limit y — > oo to get back Tgp(t), 

2B 



.a=l 



/ To(i;x) dx = -tr(SJ) . 

Jo 4 



(79) 



The significance on this term is explored thoroughly in [38]. Since it is constant it vanishes 
upon differentiation and thus makes no contribution to the vacuum energy expression. □ 
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The above method can be apphed to other integral kernels, and we can also find the 
vacuum energy density if we look at — |^T66(i; x), see [22]. 

6. Random matrix models of vacuum energy 

It has been observed by Pulling [23, 6] that level repulsion tends to decrease the magnitude of 
vacuum energy. A natural conclusion would be that in a chaotic system the vacuum energy 
should be suppressed. In this section we attempt to quantify and model this observation. 

The first serious problem is that of comparison: vacuum energy should be suppressed 
compared to what? One cannot directly compare vacuum energy of a chaotic system to that 
of an integrable one: such systems would be too different. Thus the right approach seems to 
be the average of the energy over an appropriate ensemble of chaotic/integrable systems 

In such situations it is customary to employ random matrices as models of chaotic 
systems. Which leads to a second problem: vacuum energy is not an exciting quantity when 
the spectrum is finite. Thus, (finite) random matrices do not immediately provide a suitable 
model. 

In this section we use a fusion of random matrix and graph models as a testing ground 
for the above conjecture. Namely, we study quantum graphs with equal bond lengths but 
with scattering matrices drawn from the appropriate ensembles of unitary matrices. The 
advantages are clear: each individual system will have an infinite spectrum, the spectra (in 
the limit of large graphs) will have the desired statistics and the averaging can be done 
explicitly. 

6.1. Average vacuum energy 

The spectrum of a generic quantum system with a chaotic classical counterpart is observed 
to behave like that of a random matrix, which is referred to as the Bohigas-Giannoni-Schmit 
conjecture [39]. For a system with time- reversal symmetry the appropriate ensemble of 
unitary matrices is the circular orthogonal ensemble (COE) while in the absence of time- 
reversal symmetry it is the circular unitary ensemble (CUE), which is the unitary group U (N) 
with Haar measure see [40] . For a system with time- reversal symmetry and half-integer spin 
the random matrix should be drawn from the circular symplectic ensemble (CSE). To model 
the vacuum energy of a generic chaotic system we consider quantum graphs with equal bond 
lengths where the scattering matrix S is taken from an appropriate ensemble. In drawing 
S from a random matrix ensemble we must first assume that the random matrix models a 
graph which allows scattering between all the edges, in other words a rose graph where one 
central vertex connects B loops. A random matrix S will also not, in general, correspond to 
scattering amplitudes accessible from matching conditions of a self-adjoint Laplace operator. 
However, both these limitations seem to be intrinsic to any random matrix model. Our model 
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should therefore only be regarded as a way to investigate average properties of vacuum energy 
when its generating spectrum has the desired random matrix level repulsion. 

The average vacuum energy of such a model can be evaluated using equation (22), which 
expresses the vacuum energy of a graph with equal bond lengths in terms of the eigenphases 
of S. Using the standard expression for the eigenphase density of the random matrices results 
in 



/3 



d0i... d02B = (80) 

\ z / 

j=l l<m 

where Af/3 is a normalization constant and /3 = 1 for the COE, /3 = 2 for the CUE and /3 = 4 
for CSE. 

The eigenphases of an integrable system, in contrast, behave like uniformly distributed 
random numbers on the interval [0, 2n]. This can also be modeled by integrating the vacuum 
energy expression, equation (22), over 2B independent uniform random phases, 

(^c)poisson = -Lj2^ I • • • io ^ ^^^^^Z^^) " " " ^^^^ = ^ " (§1) 

In each case the average (E^) of the vacuum energy is zero. In fact, the vacuum energy 
expression must be zero when averaged over any measure which is invariant under a rotation 
of all the eigenphases by some angle 7. Indeed, in Section 4.2 we have shown that 



1 °° 1 

^c = -^E^(t^S" + MS^n ■ (82) 

n=l 

But for any rotationally invariant distribution of eigenphases on the unit circle, (tr S") = 
for all n > and thus the mean vacuum energy is always zero. As the Casimir force is the 
derivative of Ec with respect to L this suggests that the there is no a priori reason to expect 
either an attractive or repulsive force based purely on the underlying nature of the classical 
dynamics. 

6.2. Variance 

To get a handle on how the magnitude of the vacuum energy is affected by the distribution 
of the eigenvalues we will calculate the variance of the vacuum energy for the ensembles of 
random graphs introduced previously. For Poisson distributed eigenphases the variance is 

(-B?)poisa,n = 15 I ■ ■ ■ y„ ^ BUe,/27,) d#, . . . d^B 

where 2B is the dimension of S and we used the independence of 9j to conclude that 

's2(e,/27r)S2(e,/27r)) = ^S2(^,/27r))^S2(^,/27r)) = 0. (84) 
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The variance of the vacuum energy modeled by random matrices from the circular 
ensembles can be computed using expression (30), 

(^c) = ;^ E ^;^((trS« + tr(Str)(trS™ + tr(Str)). (85) 

m,n=l 

For the circular ensembles (tr S"'tr(S''^)"^) = unless m = n as the average of a product of 
matrix elements is zero unless the number of elements of the matrix and the number from 
its Hermitian conjugate are the same [41]. Consequently 

1 °° 1 

n=l 

We notice that (| tr S"p) is the form factor of the (finite) ensemble and use the standard 



27r2L2 — 

t 

formulae [40] for CUE 




to obtain 



{2By n = 

trS"|2)QUE = \ n \n\<2B (87) 

|n| ^ 2B 

(^c)cUE = ^ {^I'^^'Km + C(3) + f *(^)(2i?)) (88) 

where ^'("^(x) is the n-th polygamma function and ( the Riemann zeta function. For all 
fixed B this is less than n'^B/SQO the variance of the Poisson distributed eigenphases. In fact 

This result parallels that for a random matrix model of the grand potential considered in 
[16]. Thus, while the variance of the Poisson ensemble grows linearly with matrix size, the 
CUE variance converges. 

The relevant parts of the form factor of the COE and CSE for finite matrix size are, 

(|trS I )cOE- AB-ny" 2B < n 

(ItrST)cSE - I '"^"^^^^^'^^ 0<n^25 

Note that in the CSE form factor the double degeneracy of the eigenphases of S (Kramers' 
degeneracy) has not been lifted. Using (86) we evaluate 

2B-1 




+ C(3)) (l + + 1/2)) - E ^^^^i^ (92) 

(93) 



2B ^(// - Z?+ 1/2) - *(„ + Z? - 1/2 
Ji^^ 2^^ 

n=2B 
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(•E')CSE = ^(Q*'"(2B)+C(3)) (l + lt{l/2-B)^ (94) 
Jf^ Mn-BjU.) ^B^,,\ (95, 

n=l / 

For S > 1 (£;e^)cOE and {E^)cSE less than (£^c)poisson- 

While the results do not have a simple closed form, in the hmit of large matrices the 
result is rather concise, 

M.JE!)cOE = ^ = tl^EDcSE (96) 

Modeling the vacuum energy variance of a quantum graph through random matrices 
suggests that the magnitude of the vacuum energy where eigenphases of S experience level 
repulsion are indeed smaller on average than those where the eiegenphases are Poisson 
distributed. Moreover, the magnitude gets smaller as the level repulsion increases from linear 
(COE) to quadratic (CUE) and quartic (CSE). Indeed, to compare the effect of the increased 
level repulsion in CSE we need to hft the Kramers' degeneracy: otherwise the degeneracy 
"compensates" the repulsion. Without the degeneracy, the result for CSE becomes 4 times 
smaller: 

Jl^^^c^^CSE/Kramers = 
thus leading to 

(^c)cOE > (^c)cUE > (^c)cSE/Kramers 

for B>1. 



7. Conclusions 



Through both the method of images and the trace formula, wc demonstrate that the vacuum 
energy in quantum graphs is a well-defined quantity (i.e. it is both convergent and a smooth 
function of the bond lengths). The closed form expression (5) is dependent only on the 
periodic paths in the quantum graph; this is a consequence of the exactness of the trace 
formula which includes only those paths. Having demonstrated how the bounce paths (closed 
paths that are not periodic) cancel when using the method of images we hope that our proof 
will shed light on the observation that periodic paths provide the correct leading asymptotic 
behavior even when the trace formula is only semiclassically correct. 

The smoothness of the expression for the vacuum energy in a quantum graph allows us 
to suggest an alternative method for its calculation, by approximating with systems with 
simpler geometries. In the case of graphs the "simpler geometry" means rational bond 
lengths, where an explicit expression for the vacuum energy is obtained. 
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We also suggest a random ensemble model for the vacuum energy when the statistics of 
the spectrum are Poisson (for integrable systems) or random matrix (for chaotic systems). 
We find the average energy in both cases to be zero, thus giving no a priori reason to expect 
a positive or negative energy from the dynamics. Furthermore, we find the variance of the 
energy and conclude that the magnitude of the energy is typically smaller when the level 
repulsion is stronger. We stress that this prediction is only correct in the probabilistic sense 
and no conclusions about particular systems can yet (if ever!) be drawn. 
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